Multi-soliton energy transport in anharmonic lattices 
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We demonstrate the existence of dynamically stable multihump solitary waves in polaron-type 
models describing interaction of envelope and lattice excitations. In comparison with the earlier 
theory of multihump optical solitons [see Phys. Rev. Lett. 83, 296 (1999)], our analysis reveals 
a novel physical mechanism for the formation of stable multihump solitary waves in nonintegrable 
multi-component nonlinear models. 



Spatially localised solutions of multi-component non- 
linear models, multi- component solitary waves, have re- 
ceived a great deal of attention in the last decade. In 
particular, recent studies in the nonlinear optics pU2] 
and Bosc-Einstein condensation [BJ have shown that, un- 
der certain conditions and only in multi-component sys- 
tems, the formation of dynamically stable localized states 
and soliton complexes is possible. Unlike their single- 
component (or scalar) counterparts, multi-component (or 
vector) solitons possess complex internal structure form- 
ing a kind of "soliton molecules" , which makes them at- 
tractive, both from the fundamental and applied point 
of view, as composite and reconfigurable carriers for a 
transport of spatially localized energy. 

Recent discovery of stable multi- component spatial soli- 
tons in optics [y|] shed a light on the general physi- 
cal mechanisms of the formation and stability of multi- 
component localised states. Such states are often called 
multihump solitons due to multiple maxima displayed in 
their intensity profile. Usually, multihump solitary waves 
appear via bifurcations of scalar solitons when a primary 
soliton plays a role of an effective waveguide ("poten- 
tial well") that traps higher-order guided modes excited 
in a complimentary field Q. On the other hand, the 
multihump solitons can be formed as multi-soliton bound 
states, when two or more different vector solitons are 
"glued" together due to balanced interaction between the 
soliton constituents M . 

Soliton bifurcations and binding enable the existence of 
multi-component localized states in many nonlinear mod- 
els; these include bound states of dark solitons Q and 
incoherent solitons |5| in optics, and multihump plasma 
waves ||. Importantly, multihump solitary waves are 
also found in higher dimensions J6|,[7|. 

The experimental and theoretical results on optical 
solitons jj],p| challenge the conventional view on multi- 
component solitary waves in other fields of nonlinear 
physics. The main question we wish to address here is: 
Can stable multihump solitons exist in other important 
models of nonlinear physics? This is a crucial issue be- 
cause, so far, the stable multihump solitons have been 



positively identified only in the nonlinear optical model 
of Refs. |l].g] that is known to possess additional symme- 
tries, which might be the reason for their unique stability. 

In this Letter, we demonstrate the existence of dynam- 
ically stable multihump solitary states in a completely dif- 
ferent (in both the physics and properties) but even more 
general model that describes the interaction of envelope 
and lattice excitations, a generalisation of the well-known 
polaron model. We reveal a novel physical mechanism for 
the formation of stable multihump solitary waves in non- 
integrable multi-component nonlinear models. 

Model. Let us consider the continuous model of the 
energy (or excess electron) transport in an anharmonic 
molecular chain |jj , described by the system of coupled 
nonlinear Schrodinger (NLS) and Boussinesq equations: 
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where t and x are the normalised time and spatial co- 
ordinate, correspondingly, ip(x,t) is the excitation wave 
function, and w(x,t) is the chain strain. The system 
is characterized by three dimensionless parameters: the 
particle mass m, the anharmonicity of the chain a, and 
the dispersion coefficient \i. 

Equation (fil) appears in a number of other physical 
contexts including, for example, the interaction of nonlin- 
ear electron-plasma and ion-acoustic waves IB], coupled 
Langmuir and ion-acoustic plasma waves [10 1, interac- 
tion of optical and acoustic modes in diatomic lattices 
p| ) particle theory models Q, etc. 

System (nl) is known to be integrable for afi — 6 p3[ . 
In this case, it possesses two types of single-soliton solu- 
tions: scalar (ijj = 0) supersonic Boussinesq (Bq) solitons 
and vector Davydov-Scott (DS) solitons Jq] which can be 
both subsonic and supersonic. Because of the complete 
integrability for afi — 6, these solitons do not interact 
with each other. For afi ^ 6, the situation changes dra- 
matically, and it has been recently shown HJlq] for the 
nearly-integrable case {afi « 6) that Bq and DS solitons 



can form a bound state for a\i > 6. From the other hand, 
it is also known that in a weakly anharmonic lattice two 
subsonic DS solitons can form a bisoliton [In]. However, 
it remains a mystery what happens when the system (TO) 
is far from its integrable limit and, especially, when the 
solitons are supersonic. In this Letter we examine the 
model (Ey) numerically for arbitrary values of afj, and, em- 
ploying the concept of soliton bifurcations, demonstrate 
the origin and exceptional robustness of multihump su- 
personic stationary solitary waves. 
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FIG. 1. Bifurcation diagram of the system (pi) for g = 7. 
Inset - close-up of the bifurcation point. Horizontal line 
- u-solitons, solid - single-hump |0, 0) solitons, dashed - 
two-hump solitons, dotted - three-hump solitons. 
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FIG. 2. Profiles of single-hump (bottom row), two- hump 
(middle row), and three- hump (top row) solitons with increas- 
ing (from left to right) values of A. 

Bifurcation analysis. The stationary solutions of Eq. 
(|l|) can be found in the form of the traveling waves 

^(x,t)^-^-e m ^ 2 - mt 4>(z), w(x,t) = -^u(z) , (2) 
4^/7 4^ 



where z = ^Jc/A[i(x — vt), and the constant c = (mv 2 — l) 
is positive for a supersonic velocity, v > 1/y/m. Substi- 
tuting Eq. (H) into Eq. (Q), we derive a system of coupled 
ordinary differential equations 

<j> zz - \<j) + 2ucj) = , 
u zz - Au + gu 2 - 4? = , (3) 

where g = ct[i is an effective anharmonicity parameter, 



and A = fi(v 2 —4fl)/c is a characteristic eigenvalue of the 
stationary localized solutions. 

Equation (g) has two types of one-soliton solutions: a 
one- component Bq soliton 



(f>o(z) = , u Q (z) = (6/5) sech 2 (z) 



(4) 



the eigenvalue problem consists of N + 
values A n = {N — n) , where n = 0, 1 



which exists for arbitrary values of g, and a two- 
component DS soliton 

4>i = 2VA(A-l)sech(VAz) , m = Asech 2 (\/Az) , (5) 

which exists only in the integrable case 5 = 6. 

To understand what happens for g / 6, we consider 
the limit (j)/u ~£« 1 and apply a multi-scale asymp- 
totic analysis. In the zeroth order in e, <j> = and Eq. 
(y) reduces to a nonlinear equation for the component 
u(z) only, with the supersonic Bq soliton solution (]4|). 
In the first order in e, we obtain a linear eigenvalue 
problem for <f>(z) characterized by the effective poten- 
tial u = uq(z). For a given value of g, the spectrum of 

1 discrete eigen- 

, ... N, and N is 

the integer part of (1/2) [^/l + (48/.g) - 1], Each cut-off 
value A„ corresponds to a bifurcation point of the node- 
less scalar soliton uq where a two-component solution 
with a nonzero component <p emerges. The latter has n 
nodes and, near the bifurcation point, can be treated as 
a fundamental (or higher-order) bound mode of an effec- 
tive potential created by the soliton uq(z). The emerging 
vector soliton can therefore be characterised by a "state 
vector" |0, n), according to the number of nodes in the 
corresponding components. 

It is easy to see that for g > 6 only bifurcations of the 
1 0, 0) state, which corresponds to the DS soliton (g), are 
possible. First bifurcation of the |0, 1) state occurs for 
the completely integrable case g = 6 at Ai = 0. In this 
case, the bifurcation pattern is identical to that of the 
completely integrable Manakov limit of two coupled NLS 
equations |Q, namely, the |0, 0) state appears at A = 1, 
and the |0, 1) state appears at Ai = 0. 

Weaker anharmonicity (smaller g) means larger num- 
ber of possible bound states supported by the effective 
potential u(z), and thus the increasing number of bi- 
furcations. Indeed, the depth of the effective trapping 
potential is inversely proportional to g. The |0,0) state 
always exists, even for a shallow potential u(z). 

We now consider in detail the formation of multihump 
solitons in the cases of weak {subcritical, g < 6) and 
strong (supercritical, g > 6) anharmonicity, respectively. 

Supercritical regime. In the absence of bifurcat- 
ing higher-order solutions, the multihump solitons are 
formed only via binding of the |0, 0) vector solitons. The 
physics of this mechanism is simple. The interaction 
forces between closely separated fundamental solitons are 
different for both the 4> and the u components. Namely, 



while two in-phase u-solitons attract, the two in-phase <f>- 
solitons repel. This allows for the existence of multihump 
nodeless modes of the field <\> trapped in the multi-well 
potential u. Each of such multihump solitons can be con- 
sidered as a bound state of several |0, 0) DS solitons, with 
in-phase humps in both components. 

It is convenient to represent the solution families as 
branches on the bifurcation diagram Q vs. A, where 
Q = Q u + Q,p = J u 2 dz + J (jPdz is the total soliton 
power. Typical bifurcation diagram for a supercritical 
case (g = 7) is shown in Fig. RL Solid line represents 
the bifurcating solution |0,0), and it can be seen on the 
close-up of the bifurcation region, that the branches rep- 
resenting two- and three-hump solutions start off at the 
same point A = Ao but with the energies approximately 
equal to that of two or three u-solitons. Examples of such 
multihump solitons are shown in Fig. [2|, and it is clear 
that this novel type of soliton bifurcations occurs from a 
countable set of infinitely separated single solitons. With 
increasing A, separation between the humps decreases un- 
til all solitons of this type become single-humped (Fig. @, 
right column). 
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FIG. 3. Bifurcation diagram of the system (pi) for g — 1. 
Shown are branches of the fundamental u— soliton (hor- 
izontal), |0, 1) (dot-dashed), |0,3) (dotted), |0, 2} (solid), 
and |0,4) (dashed) solitons. Below: profiles of |0, 1} and 
|0, 2) solitons (top row) and some examples of their bound 
states, |0,3) = |0, 1) + |0, 1) + |0, 1) and |0, 4} = |0, 2) + |0, 2} (bot- 
tom row). 

Subcritical regime. In this case, bifurcations of the it- 
state do not lead to multihump solitons. That is, in sharp 
contrast to the coupled NLS equations describing vec- 
tor solitons in nonlinear optics, none of the higher-order 
states |0,n) become multihumped in the system under 
consideration. Although the function cf> does have mul- 
tiple maxima in its intensity profile, because of the non- 
self-consistent source for the u-component, it does not 
cause significant distortions in the shape of the effective 
potential, u(z), and the total intensity I{z) = u 2 + (j) 2 re- 



mains single-humped. Typical bifurcation diagram for 
the case g = 1 is presented in Fig. g. In this case 
N = 2, and only bifurcations of the |0, 1) (dot-dashed 
line) and |0, 2) (solid line) solitons are shown. Corre- 
sponding modal profiles of the bifurcating solitons are 
presented in Fig. || (top row). 

Similar to the supercritical regime, the multihump soli- 
tons can exist only as bound states of the bifurcating |0, 1) 
or |0,2) solitons. They appear at the bifurcation points 
A n , and they have energies equal to a number of lower- 
order solitons "glued" together by the low-amplitude 
components. The number of nodes that the </>-component 
has in the composite soliton depends on the number of 
|0,n) solitons forming that bound state. Typical exam- 
ples of such solutions are presented in Fig. (bottom 
row). 

From this analysis, we can conclude that, in this model, 
the multihump solitons appear as bound states of |0, 0) 
solitons for any value of anharmonicity parameter g. 
In addition, multihump solitons of more sophisticated 
modal structure are also possible. 
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FIG. 4. (a,b) Stable dynamics of a three-hump soliton for 
g = 7 and A = 0.824. The family of these solitons is shown 
in Fig. hi (dotted line), and the corresponding profiles are 
presented in Fig. H (top row). (c,d) Stable dynamics of a 
four-hump soliton for g = 3.005 and A = 2.9. 

Dynamical stability. The second equation of the sys- 
tem (l|) is the so-called "ill-posed" (or "bad") Bq equa- 
tion [18[. It possesses an intrinsic linear instability and 
therefore its reliable numerical solution for non-zero /i 
is unfeasible. This linear instability is not inherent in 
the original physical model, and it can be traced to ne- 
glecting higher-order spatial derivatives in Eq. (pi). In 
the case of the energy transport in anharmonic molecu- 
lar chains, Eq. (Q) with /i = 1/12 originates from the 



following system of discrete equations 1 14 



11 



i^ n = -A 2 (tf n ) - (W n + W„_i)tt„ , 

at 



:(l* 



n+l\ 



I*. 



I*, 



I* 



n+2 



(6) 



I 2 ), 



where A.2(X n ) = X n+ i+X n _i — 2X n . Discrete functions 
ty n and W n define, in the continuous limit, the excitation 
wave function, ip(x,t), and the strain function of the lat- 
tice, w(x, t). 

Therefore, it would be justified to study the dynamics 
of the stationary solutions of Eqs. (||) numerically by em- 
ploying the original discrete dynamical system (ra). Be- 
sides, the argument can be reversed, and such a discrete 
system can be treated as a regularised numerical discreti- 
sation scheme for a general system of coupled NLS and 
ill-posed Bq equations. 

We investigate the dynamical stability of the mul- 
tihump solitons for two distinct cases of a subcritical 
and supercritical anharmonicity discussed above. The 
condition of a unit norm for the envelope function ^> n 
(S n l*n| 2 = 1) i s satisfied in all cases, and v is chosen 
close to the sound velocity (c = 0.0244) to allow for a 
smooth discretisation. 

In the supercritical regime (g > 6), where multihump 
solitary waves can be formed through binding of several 
|0, 0) states together, our numerical simulations indicate 
that such solitons are stable as long as the separation 
between the humps is sufficiently large. This property is 
just opposite to that observed for multihump optical soli- 
tary waves |2| . An example of the stable dynamics of a 
three-hump soliton for g = 7 is shown in Figs. 0(a,b). 
All solitons of the DS type, i.e. |0, 0) soliton states, pre- 
sented in Fig. n] (by solid line) and Fig. (bottom row), 
exhibit similar stable dynamics. 

It is important that the same mechanism of the cre- 
ation of multihump solitons applies to the subcritical 
regime. This means that the dynamically stable multi- 
soliton bound states described above exist also for g < 6. 
As an example, propagation of a four-hump soliton at 
g = 3.005 is demonstrated in Figs. [|(c,d). After initial 
adjusting of the soliton amplitudes (due to the discreti- 
sation), only small amplitude breathing occurs [see Fig. 
f|(d)], otherwise the soliton dynamics is stable. In con- 
trast, all bifurcating higher-order solitons are dynamically 
unstable. 

Our results on the robustness and stability of multi- 
soliton states call for a systematic revision of our un- 
derstanding of the role of nonlinear localized modes in 
a number of physical phenomena related to the nonlin- 
ear transport in macromolecules B and even artificial 
nanoscale structures p9f , where the coupling of two (or 
more) degrees of freedom occurs. How the soliton binding 
and existence of multi-soliton states modifies the nonlin- 
ear kinetics, nonequilibrium thermodynamics 2C }, and 



strated their dynamical stability. Along with the recent 
studies on multihump optical solitons [0J|], these results 
call for re-examination of the role of multi-component 
solitary waves in other fields of nonlinear physics. 
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Fund of The Australian National University. 



other properties of the system! These questions remain 
to be answered. 

In conclusion, we have found robust two-component 
solitary waves in a polaron-type model of the energy 
transport in anharmonic lattices. We have revealed a 
novel physical mechanism for the formation of multihump 
solitons in a discrete anharmonic lattice and demon- 
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